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We investigate the probability distribution of the quantum fluctuations of thermodynamic func- 
tions of finite, ballistic, phase-coherent Fermi gases. Depending on the chaotic or integrable nature of 
the underlying classical dynamics, on the thermodynamic function considered, and on temperature, 
we find that the probability distributions are dominated either (i) by the local fluctuations of the 
single-particle spectrum on the scale of the mean level spacing, or (ii) by the long-range modulations 
of that spectrum produced by the short periodic orbits. In case (i) the probability distributions are 
computed using the appropriate local universality class, uncorrelated levels for integrable systems 
and random matrix theory for chaotic ones. In case (ii) all the moments of the distributions can 
be explicitly computed in terms of periodic orbit theory, and are system-dependent, non-universal, 
functions. The dependence on temperature and number of particles of the fluctuations is explicitly 
computed in all cases, and the different relevant energy scales are displayed. 

I. INTRODUCTION 

Self-consistent theories which admit almost independent motion of quasiparticles and give rise to Fermi gases 
constitute a basic tool to describe fermionic many-body systems. Well known examples are independent-particle 
models in nuclear physics, quantum dots and wires in condensed matter physics and, within a mean- field theory for 
the valence electrons, the jellium model for the electronic properties of metals and of simple alkali metal clusters. 
One of the most spectacular predictions of mean-field theories is the occurrence of shell effects in the energy of the 
fermion gas as a function of the number of particles, that leads to magic numbers and deformations in nuclei and 
metallic particles The shell and super-shell structures are quantum fluctuations that describe the bunching of 

the single-particle energy levels with respect to their average behavior. From a scmiclassical point of view 1^,^, the 
bunching of the energy levels is interpreted as a modulation produced by the classical periodic orbits. 

Similar quantum fluctuations are present in all the equilibrium thermodynamic functions describing the many-body 
system. As we will see in detail in the following, there are two main features of the single-particle spectrum that 
influence the thermodynamic quantum fluctuations of a Fermi gas. On the one hand, the spectral fluctuations on local 
scales, i.e. in the range of the mean level spacing 6 between single-particle energy levels. These local fluctuations are 
universal, and the universality class depends on the nature of the single-particle dynamics ||,^. There are two main 
classes, corresponding to regular and irregular (fully chaotic) motion, respectively (with subclasses in the latter case 
depending on the presence or not of time reversal symmetry). The local statistics of regular systems are described by 
a set of uncorrelated levels, while the chaotic ones by random matrix ensembles. 

The second ingredient that strongly influences the thermodynamic fluctuations are the shell effects, i.e. the long 
range modulations of the spectrum produced by the short periodic orbits. These occur on a scale Ec related to the 
inverse time of flight across the system, with typically Ec ^ S. Contrary to the local fluctuations, the long-range 
modulations are system-dependent, since the short orbits are. The amplitude of this effect depends on the nature of 
the dynamics. If the classical motion is regular, symmetries are present, periodic orbits come in degenerate families, 
and shell effects are important. On the contrary, if the corresponding classical motion is chaotic, the periodic orbits 
are isolated and the long-range modulations will have, compared to the regular case, a smaller amplitude. Different 
aspects of both quantum phenomena, namely the universality of the single-particle fluctuations on a scale d and 
the appearance of long-range shell effects on scales Ec, and their influence in the thermodynamics of Fermi gases, 
have been considered in nuclear and atomic physics, and in different condensed-matter systems (see for example 
Refs. [0j|J^,P). In most cases the interest concentrated on particular systems, or in the study of mean properties. 

Our main goal is, within a Fermi gas approximation, to provide a statistical analysis of the quantum fluctuations 
of the different thermodynamic functions, like for example of the total energy of the gas, the entropy, etc, and to 
explicitly compute their probability distribution function. Special attention is devoted to the temperature dependence 
and to the typical magnitude of the fluctuations. Other aspects, like the autocorrelation functions and the connexion 
between the fluctuations of the gas and those of the single-particle levels at Fermi energy, are also investigated. We 
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concentrate in the fluctuations of equilibrium properties of ballistic Fermi systems (i.e. with no disorder present in 
the potential) whose mean-field single particle motion is either regular or fully chaotic. Schematically we have in mind 
a mean-field potential which is almost flat in the interior and rises at the surface, like the Woods-Saxon potential in 
nuclear physics, or the static homogeneous potential produced by the ionic background in the jellium model. Under 
these conditions, and in particular in the extreme case of a sharp boundary, the nature of the classical dynamics is 
determined by the shape of the cavity and by the external fields acting on the system. We will ignore physical effects 
related to the spin degree of freedom, and we will moreover assume that the inelastic processes can be ignored, namely 
the coherence length of the fermions is much larger than the system size. 

The picture that emerges from our analysis is very rich. The quantum thermodynamic fluctuations can be either 
dominated by the local universal or by the long-range system-specific fluctuations of the single particle spectrum. 
Among these two possibilities, the prevailing aspect depends on the quantity considered, on the nature of the un- 
derlying classical dynamics, and on temperature. When the fluctuations are dominated by the local fluctuations, 
the corresponding thermodynamic probability distributions are universal. When conveniently normalized, they only 
depend on the temperature and the quantity considered, but not on any speciflc property of the system (aside the 
regular or chaotic nature of the single-particle motion). The universality of the thermodynamic fluctuations follows 
the classification of the local fluctuations given above: one universality class for systems whose classical single-particle 
dynamics is regular, another for those with a chaotic dynamics (with subclasses depending on the presence or not of 
time reversal symmetry). Universal fluctuations are observed, for some thermodynamic functions, at temperatures 
much smaller than Ef,. Examples are the fluctuations of the number of particles for systems with flxed chemical 
potential when the underlying classical dynamics is chaotic, and those of the entropy and the speciflc heat for reg- 
ular and chaotic motion. The corresponding probability distributions are found to be, in general, non-Gaussian. 
For temperatures of order 6 their shape is very sensitive to temperature variations, but they always remain in the 
corresponding universality class. When temperatures of order Ec are reached, the universality is lost. At those tem- 
peratures the short orbits become dominant, producing system-dependent effects. The moments of the distributions 
can then be computed from periodic orbit theory, as described below. At still higher temperatures the fluctuations 
are exponentially suppressed. 

In contrast to the flrst type described above, sometimes the thermodynamic fluctuations arc dominated by the long 
range modulations of the single-particle spectrum produced by the short periodic orbits. Due to the system speciflc 
nature of the short orbits, the corresponding probability distributions are non-universal, non-Gaussian functions. 
We have obtained explicit convergent expressions to compute all the moments for this type of fluctuations in terms 
of periodic orbits. Interference effects between repetitions of the orbits produce, in particular, non-vanishing odd 
moments and lead to asymmetric distributions. The fluctuations of all the thermodynamic functions are of this type 
for temperatures of order Ec or higher. However, some thermodynamic functions have non-universal fluctuations at 
any temperature. Examples of the latter behavior are the fluctuations of the total energy of the gas, the response 
of the system when an external parameter is varied and, only in the case of a regular single-particle dynamics, the 
fluctuations of the number of particles for systems with fixed chemical potential. Due to the dominance of the long 
range modulations, the probability distributions are relatively insensitive to temperature variations up to temperatures 
of order Ec- At temperatures higher than Ec the typical size of the fluctuations decay exponentially. Finally, the 
autocorrelation functions in this case present long range order. 

As a general rule, but with exceptions depending on the quantity and on the temperature considered, clear signatures 
of the regular or chaotic nature of the dynamics are found in the fluctuations. For example, their typical size is in 
general much larger in integrable systems compared to chaotic ones. The relative amplifying factor is found to be some 
growing function of the adimensional parameter g = Ec/S. Besides, and irrespective of the nature of the dynamics, 
the fluctuations of the gas are always much larger than the corresponding fluctuations of the single-particle energy 
levels at Fermi energy. Their relation is also expressed in terms of g. 

The paper is organized as follows. Section || contains the basic definitions concerning the thermodynamics of a 
Fermi gas. An arbitrary thermodynamic function is decomposed into a smooth plus an oscillatory part. The smooth 
part describes the usual thermodynamics of the gas (bulk properties), while the oscillatory part are the quantum 
finite-size fluctuations. The latter are expressed in terms of the periodic orbits of the classical motion associated with 
the corresponding single-particle dynamics. The validity of this approximation and the main energy scales involved in 
the problem are also introduced. In Section [II the variance of the thermodynamic functions is written in terms of the 
spectral form factor, i.e. the Fourier transform of the single-particle spectral two-point correlation function (Eq.(p9|)). 
This connexion allows to make a qualitative analysis of the behavior of the moments of the probability distribution 
of the different functions and for the different kinds of classical motions. The second moment is explicitly computed 
in Section IV using a schematic approximation for the short-time dynamics that provides simple but qualitatively 
useful expressions. The dependence on particle number and on temperature is explicitly considered. The complete 
characterization of the probability distribution of the fluctuations is included in Section The moments are explicitly 
computed for quantities whose distribution is dominated, by the short orbits. On the contrary, local uncorrelated levels 
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and random matrix fluctuations are used to compute the distribution of universal functions (at low temperatures). 



Similar arguments are exploited in Section Vl to compute the autocorrelation of the thermodynamic functions as 



a function of the different parameters. The application of our general results to d iffer ent systems, like the atomic 



nuclei, mesoscopic quantum dots and metallic clusters is briefly discussed in Section VII. This section contains also a 



discussion on the relation between the fluct uatio ns of the gas and those of the individual single-particle energy levels. 



Our final remarks are included in Section VIII. All over the paper we compare some of the results with numerical 



simulations of Fermi gases contained in regular and chaotic cavities. 

II. GENERAL SETTING 
A. Basic equations 

When the chemical potential /x is fixed the thermodynamic behavior of a non-interacting Fermi gas is described by 
the grand potential 

n{x, M, T) = -k^T J dE p{E, x) log[l + e^P-E)/kBT^^ ^ (1) 

where 

p{E,x)^g,Y,6[E-E,{x)] (2) 

3 

is the density of states of the single-particle energy levels Ej{x). The parameter x denotes the dependence on some 
external parameter, T is the temperature, and is Boltzmann's constant. The prefactor = 2 when spin degeneracy 
exists. For simplicity we will from now on omit it, but it can be easily restored in the final results. 
When the number of particles N in the gas is fixed, the energy of the system is 

U{x, N,T)^ JdEE p{E, x) fiE, ^, T) , (3) 

where fi is determined from the equation 

N = J dE p{E,x) f{E,^Ji,T) , (4) 

and / is the Fermi function 

Knowing these two thermodynamic potentials, other thermodynamic functions are directly obtained by differenti- 
ation. For example, in the grand canonical formalism the expectation value of the number of particles in the system 
is 



(6) 



Another important thermodynamic quantity is the reaction of the system to a variation of the external parameter a;, 
characterized by the response function (here in the grand canonical formalism) 



R{x,p,T) = - ^ 



(7) 



The physical interpretation of R depends on the sample geometry and on the nature of the parameter x. R is the 
orbital magnetization when x is a magnetic field, a persistent current when the sample geometry is annular and x is an 
Aharonov-Bohm fiux, or a force (pressure) when x controls the shape of the confining potential (like when deforming 
a nucleus for example; at equilibrium, the nucleus shape is defined by i? = 0). Finally, we will also consider the 
entropy 
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df 



(8) 



from which, in particular, the many-body level density may be computed by exponentiation, and the specific heat (in 
the canonical ensemble) 



dU 
dT 



(9) 



The case of the susceptibility, i.e. the variation of the response function (the second derivative of the grand potential) 

(10) 



x(^,M,r)= — 



is a more subtle problem and will not be considered here in detail. 

Assuming the single-particle spectrum Eq. (^) known for any x, the integration in Eq.([^) is straightforward. It leads, 
by differentiation, to sums over the single-particle energy levels whose structure depends on the quantity considered. 
As the temperature is lowered the Fermi factor f{E,iJ,,T) involved in the sums tends to the Heaviside step function, 
which restricts the sums to the single-particle energy levels satisfying the condition Ej(x) < fi. In this limit the grand 
potential is closely related to the ground state energy of the Fermi gas (it is the sum of the single particle energy levels 
measured with respect to the Fermi energy), J\f is the number of occupied single-particle levels, and the response 
of the system is the sum up to the Fermi energy of the single-particle contributions dEj/dx. They all involve the 
ensemble of particles of the Fermi gas. In contrast, other quantities, like the entropy and the specific heat, vanish at 
zero temperature. 



B. Semiclassical approximation 



We now express the thermodynamic functions in terms of classical quantities related to the single-particle dynamics 
(cf Ref. ^ for a good introduction). The density of states in Eq.(|l|) is written as a sum of smooth plus oscillatory 
terms 

p{E,x) ^p{E,x) +p{E,x) . (11) 

This approximation is valid in the semiclassical regime 5 where the typical wavelength of the fermions at energy 
jjL is much smaller than the system size, and fails at the bottom of the spectrum. The first, smooth term p corresponds 
to the bulk contribution. For example, when the gas is simply contained in a finite volume V of space having an 
arbitrary shape, and the single-particle levels are filled up to a Fermi energy Ep with particles of mass m, the average 
number of occupied levels is, to leading order 

AT^i^f^V'^r, (12) 
while the mean level density p and corresponding mean level spacing 5 between single-particle states are 



2V 



3/2 



eI'^ . (13) 



The substitution in Eqs.(0) and (||) of the density of states by its average behavior, Eq.(p^, gives rise to the smooth, 
bulk classical expressions for the thermodynamics of the gas [llo| . In finite fermionic systems, like the atomic nucleus 
or a metallic particle, the smooth part coming from the single-particle potential is of no use, and is rather computed 
in a self-consistent manner from the original many-body problem [ prf . 

The most relevant contribution to the many-body problem is the fluctuating part, which is of quantum mechanical 
origin, and describes the discreteness of the single-particle spectrum. It is given by a sum over the classical periodic 
orbits |^,|l^. To leading order in an ?i-expansion 

oc 

p{E,x) ^2^^Ap,r{E,x)coa[rSp{E,x)/h + iyp,r] . (14) 

p r—1 



4 



The sum is over all the primitive periodic orbits p (and its repetitions r) . The orbits are characterized by their action 
Sp, period Tp = dSp/dE, stabiUty ampUtude Ap,r, and Maslov index Vp^r- 

ExpUcit expressions for Ap^r are available. Their functional form depends on the nature of the dynamics. For a 
d-dimensional integrable system, where action-angle variables exist, the periodic orbits are organized on resonant 
tori, i.e. d-dimensional manifolds having the topology of a torus with commensurate frequencies. If Ij are the action 
variables, then the frequencies are Uj = dH/dlj = 27rmj/Tp, where the rrij are integers that label the periodic orbits 
(that we abbreviate in the unique index p). When the {mj} have no common divisor they label a primitive periodic 
orbit, otherwise it is a repetition. The Ap's take in this case the form H 

Al = . (15) 



p 



det{du;j/dlk}p J2j ' dlj/dr,, 



The repetitions do not appear explicitly, and are considered as primitive periodic orbits. Instead, in chaotic system 
the amplitudes are (l2j 

Ap,r = , ""^ = . (16) 

IdetfAfr- ni 



Here / is the identity matrix and Mp the monodromy matrix obtained from the linearization of the equations of 
motion in the vicinity of the corresponding primitive periodic orbit. In both cases the Maslov indices i^p ^ are constant 
phase factors (locally) independent of the energy and of the external parameter. 

Eq.(^^ is inserted in Eq.(|^) to compute the oscillatory part of the grand potential. To leading order in h and for 
low temperatures {kgT <C n, degenerate gas approximation) 

h{x, n, T) ^2f? ^"''^^^l'" """^ ^os{rSp/n + Vp.,) . (17) 

p r—l P 

The classical functions entering in this expressions are evaluated ai E = fj,, and they all depend on the external 
parameter x. The temperature introduces the prefactor Kt{t), which acts as an exponential cut off for the long orbits 



'^-W = • w^T , , = h/{2TT^k^T) . (18) 

smh(T/rT) 

For temperatures such that Tr^ <C Tmi„, with Tj^in the period of the shortest periodic orbit, the quantum fluctuations 
are washed out and only the smooth behavior given by p survives. 

Similarly we may analyze the energy of the gas at a fixed number of particles. The smooth part of U differs from 
that of n, but to leading order in a semiclassical expansion we find that their fluctuating parts coincide. U is thus 
given by an equation identical to Eq. with /i not a continuous variable but now a function of N, x and T 

determined by Eq.(^ but where p is used instead of p, 

U{x, N, T) = Q.{x, p{x, N, T),T) . (19) 



The prefactor Ap^r k,t{t Tp)/Tp in Eq. ( |17| ) varies smoothly with the external parameter, the chemical potential, 
or temperature. The main contribution to the oscillatory behavior of the grand potential comes in fact from the 
variations of the phase factor rSp{x, p)/h with respect to x and p. When differentiating Eq. ( p7| ) or Eq. ( p^ ) to 
compute the fluctuating part of other thermodynamic quantities according to Eqs. (|6|)-(p^ we only keep, to first order 
in a semiclassical expansion in h, the terms coming from the variations of the phase factor. In this approximation, 
the fiuctuating part of all the thermodynamic quantities considered has the same structure, and can be written in a 
compact form as 

oo 

^{x,p,T) ^2 C^^Ap^r{x,p,T)cos[rSp{x,p)/n + iyp^r] ■ (20) 

p r—l 

The constant coefficient C includes all the terms not depending on the periodic orbits. This coefficient, as well as 
the orbit-dependent prefactor Ap^r computed in the leading semiclassical approximation, are given for some of the 
thermodynamic functions in Table I. When necessary, the sine functions have been systematically transform into 
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cosines by subtracting a 7r/2 to each Maslov index. The prime in n'^ denotes the derivative of with respect to 
temperature, k'^ = dKr/dT, while 



dx 



(21) 



Eq.(gO|), together with Table I, are the basic expressions on which the study of the probability distributions of the 
fluctuations is based on. The knowledge of the periodic orbits (or at least of those satisfying, at finite temperature, 
Tp ~ Tj-) allows to explicitly compute, using Eg. (pO|) , the quantum oscillatory contributions to <i>. We are interested, 
however, in a statistical analysis valid for a generic system rather than in an explicit computation for a particular 
one. 





C 


•^p.r 


F{t) 


u 




Ap^r l^rij'Tp) / r^T^ 


4(r) 


N 


n 


Apr KrpirTp) / TTp 


t2 4(r) 


R 


n 


Ap r Qp K-j-irTp) / TTp 




§ 




Ap.r (^''p)/^^'''p 





Table I: the coefficient C and prefactor Ap^r of Eq. (po|) for different thermodynamic functions. Sec the text for details 



In the third column is listed the function F(t) appearing in Eq. (29) 



As T ^ the sum (po| ) over periodic orbits is in general divergent. For instance, in fully chaotic systems it diverges 
for temperatures T ~ h^/2nkB (where ^ is the Lyapounov exponent), and converges otherwise. Using resummation 
techniques finite sums may be obtained, as was done for example for the magnetic susceptibility of two-dimensional 
billiards in However, we will not need to proceed along these lines since, as we will see, the moments of the 

thermodynamic quantities based on these sums are, in some cases, convergent. In the cases where they don't, the 
distributions will be directly determined from the local universal fluctuations. 



C. Averages, energy and parameter scales 



There are several different energy scales that are relevant when considering the quantum thermodynamic fluctua- 
tions. The smallest energy scale in the problem is the single-particle mean level spacing S ^ p^^ at the Fermi energy 
fi (since fc^T <C /x, we don't make in this qualitative considerations the difference between the chemical potential and 
the Fermi energy). The time scale associated to this energy is the Heisenberg time 

T„ = h/S . (22) 

The largest energy scale is the Fermi energy Ej,, or chemical potential /i. The third relevant energy scale has a simple 



semiclassical origin. According to Eq. (20), locally each periodic orbit p produces a periodic modulation (in energy) 
of wavelength Xp ^ h/Tp{fi). Since there are orbits of arbitrary long period, Xp can be arbitrarily small. There is on 
the contrary an upper bound to Xp determined by the period Tmin of the shortest periodic orbit (at E = p). The 
corresponding energy scale is 

Ec = h/Tmin ■ (23) 

This is the largest scale in which long-range modulations of the single-particle spectrum occur, and flxes the size of 
the shell structures at E = Ep. The ratio of E^ to S is another important parameter 

g = r„/T™,„ Ec/5. (24) 
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The energy and the adimensional parameter g are, in some respects, the analog for ballistic systems of the Thouless 
energy and the adimensional conductance defined in diffusive systems, respectively Its physical meaning in 

ballistic systems is clear from Eq. (|^. It is the number of single-particle states contained in the last shell, which is 
of size Ec- It is therefore a measure of the collective effect of the modulations produced by the shell contributions, a 
"shell-strength" parameter. 

Assuming the typical system-size is L, it is useful to display the different energy scales 5, Ec, and /i in terms of the 
typical number of the De Broglie wavelengths the system can accommodate at Fermi energy, kpL, where kp is the 



Fermi wavevector. These estimates are based on the generalization of Eqs. and (13) to a d-dimensional cavity. 
On the one hand we have g — Ed 5 cx {kpLY~^. Moreover, ^/Ec — kpL/2TT. Thus, in the semiclassical regime 
kpL 3> 1 where our description applies, the parameter g is large, and the different energy scales are well separated, 
S <^Ec<^ fi. 

The fourth relevant energy scale is the temperature. Eq. (|2^) is valid in the limit kgT ^ /i of a strongly degenerate 
gas, which is the relevant limit for most physical applications. Notice the large coefhcient 2tt'^ « 20 involved in the 
parameter in Eq.(|l8|), as compared to Eqs. ( p^ ) and (p^). This large coefficient cannot be ignored, and therefore 
the relevant thermal energy to be compared with the quantum scales Ec and S is 2iT^kBT . 

The fluctuating part of a given thermodynamic function $ shows, as a function of the chemical potential ^, 
oscillations described by Eq. (pO|). The statistical properties of and in particular its probability distribution, will 
be computed in a given interval of size A/x around /i. This interval must satisfy two conditions. It must be sufficiently 
small in order that all the classical properties of the system remain almost constant. This is fulfilled if A/i <C /i. 
Moreover, it must contain a sufhciently large number of oscillations to guarantee the convergence of the statistics. As 
stated previously the largest scale associated to the oscillations is E^- Then clearly we must have A/i ^ E^,. In the 
semiclassical regime the hierarchical ordering between the different scales is therefore 

(5 < < A/i < /i . 

Since ji/Ec oc kpL, a typical scale for the smoothing energy window is Afi/Ec oc {kpLY^^. We thus define the energy 
average of a certain oscillating function as 

(/(m))^ ^ / fWfi' . (25) 



In a similar way a parameter average is defined as 



1 



x+Ax/2 



{f{x))x ^ / f{x')dx' , (26) 

L^X Jx~.Ax/2 

where Aa; is defined by similar arguments as A/i. As in the energy average, the classical properties must not signif- 
icantly change in the interval Ao;, but it must contain several oscillations to make a statistical analysis appropriate. 
The largest scale associated to these parameter oscillations is again related to the shortest orbit, and given by h/Qmin- 
We finally note that if all these conditions are satisfied both the energy and parameter averages gives the same 
result, and consequently commute. 



III. QUALITATIVE ANALYSIS OF THE DISTRIBUTIONS 

By definition, the average value of the fluctuating part $ is zero. Aside the average value, the variance is the more 
basic aspect of the probability distribution of the fluctuations. It provides the typical size of the oscillations, and can 
easily be compared with experiments. We will now compute a general expression for the second moment that allows 
to make a qualitative analysis of the distribution of the different thermodynamic quantities. Similar arguments for 
the higher moments lead to similar conclusions. 

From Eq. (^o|) the square of $ is expressed as a double sum over the periodic orbits involving the product of two 
cosines. The latter may be expressed as one half the sum of the cosine of the sum and that of the difference of the 
actions. The average over the term containing the sum of the actions vanishes, due to its rapid oscillations on a scale 
A/i. Therefore 



<^^)^2C'{J2A,A,,cos(^^^)) . (27) 
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To simplify the notation we have momentarily considered the repetitions of a primitive orbit as a different primitive 
periodic orbit, and have included the Maslov indices in the definition of the action. Ordering the orbits by their period, 
and taking into account the restrictions imposed by the averaging procedure, we can relate the variance Eq. ( p7| ) to 
the semiclassical definition of the form factor K{t) (i.e., the Fourier transform of the two-point correlation function 
{p{^ + e/2)p{^ — e/2))^ with respect to e), expressed as 



K{t) 



^ Ap Api cos 



Sri 



(Tp + Tp,) 



(28) 



The form factor defined in Eq. (28) has units of time. When expressed in terms of K{t), the variance of the quantum 
thermodynamic fluctuations considered above takes the simple form 



($2) = 



C2 



271^ 



dT 



K{t) F{t) , 



(29) 



where F{t) depends on the thermodynamic function considered (cf Table I). Analogous expressions connecting the 
variance of different quantities like the persistent current or the conductance fluctuations have been obtained when 
describing disordered metals [0. In the latter case, the form factor for times r <C t„ is closely related to the classical 
return probability for a diffusive particle. 



To obtain Eq.(|29|) we made use of the fact that the orbits giving a non-zero contribution to (27) have similar actions 
(unless their average will be zero). This implies that their period is also similar, and can be considered to be the 
same in the prefactor (but not in the argument of the oscillating part). Similar arguments apply to the Qp's. It 
can be shown by arguments invoking the stationarity of the form factor under variations of the external parameter x 
that the contributing orbits have very similar derivatives Qp. The factor Q^(t) entering the variance of the response 
function has a statistical meaning. When replacing the double sum in Eq. (|2^) by an integral over the period t, among 
all the periodic orbits p of period between t and t + dr we consider the distribution of the derivatives (^). Then 
Q^{t) is the second moment of that distribution. The expression ( p9[ ) for the response function therefore neglects the 
non-statistical behavior of the short orbits. 

By definition, the form factor ( p8| ) is system dependent. But general statistical statements, depending only on the 
nature of the underlying classical dynamics, have been conjectured for ballistic systems. When the classical dynamics 
is integrable, the quantum energy levels are believed to behave as an uncorrelated sequence This implies no 
r-dependence of the form factor. 



Kint (t) = Tfl 



(30) 



where t„ was defined in Eq. (|22D. The situation is different for fully chaotic systems, where it has been conjectured 
that the fluctuations are described by the corresponding ensembles of random matrix theory For the form 

factor it gives 



Krmt{T) 



2r-Tlog(l + |j) 



T e(T„ -r) +r« e(r-T«) 



2th - T log 



2r 



2r 



(3 = 1 
/3 = 2 



(31) 



where /? = 1 (2) for systems with (without) time-reversal symmetry (since we are neglecting the role of spin we do 
not include the symplectic symmetry). For short times t <ti Th this function behaves as 



(32) 



Aside the form factor, to compute the variance of the response function information about the Qp's is required. 
For chaotic systems numerical evidence as well as general considerations suggest a Gaussian distribution, with zero 
average (if the area of the cavity is preserved when x varies) and variance proportional to the period t . Then 
Q'^{t) ~ ar, where a is some system-dependent constant which corresponds to a classical diffusion coefficient pO[ |. In 
the integrable case no universal distribution exists for the Qp's, but explicit calculations for some integrable systems 
as well as general heuristic arguments indicate that the second moment is proportional to [ pT| . In summary 



ar chaotic 



lut'^ integrable . 



(33) 
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For short times, Eqs. (|3^) and ( |3^ ) are good approximations in real systems in the regime Tmin t <^ Th- Their 
semiclassical origin is related to the statistical behavior of the long classical periodic orbits 22 1^. 

For times of the order of Tmin no statistical universal behavior of the form factor exists, and Eq. (|3^) and ( |3^ ) 
do not provide a good description. For such times the off-diagonal contributions in Eq. ( P8|) are eliminated by the 
averaging procedure, and the behavior of K{t) is described by a series of delta peaks at r = Tp obtained from the 
diagonal terms p — p' in the double sum 

K{T)^h'Y.KHr-rj,) . (34) 
p 

The lowest peak is located at t = Tmm, and for r < T„iin the form factor is identically zero. 

Based in Eqs.(p9|)-(^4|) we now make a simple qualitative analysis of the variance of the fluctuations of the thermo- 
dynamic functions. To simplify the analysis, in the chaotic case we will use the simple approximation ( ^2|) extended 
to times r —^ cjo, and ignore the saturation of the form factor for longer times described by Eqs. (^l]). By doing this, 
we are overestimating the weight of the long orbits. To begin with, we also set the temperature to zero. In Table I 
we put Kt — 1 and k'j, = 0. The variance of the entropy fluctuations therefore vanishes, as it should. We now discuss 
the non-vanishing quantities. 

Consider for example the variance of the energy (or of the grand potential). From Eqs. ( |3^ ) and ( p3| ) we see that for 
chaotic systems the integrand in Eq. (^oj) behaves as t~^, while from Eqs. (po|) and (^) it follows that for integrable 



motion the integrand varies as t~*. Therefore in both cases the integral (29) converges for long times, even at zero 
temperature. The dominant contributions come in fact from short times, where the divergence of the integral is 
stopped by the cutoff at r = Tmin introduced by the form factor. When applied to other thermodynamic quantities, 
the same analysis shows that the fluctuations of the response function for integrable and chaotic dynamics, and those 
of the particle number for integrable systems, are also dominated by the short-time contributions. Because the short- 
time dynamics is specific to each system, then in general the second moments of those functions are non-universal, 
and consequently the same is true for the probability distribution. This conclusion remains valid as the temperature 
increases, since a finite temperature truncates the long orbits. 

We thus conclude that the Ructuations of the grand potential, the energy and the response function for integrable 
and chaotic dynamics, and those of the particle number for integrable systems, are dominated, at an arbitrary 
temperature, by the shortest non-universal periodic orbits of the system. 

As mentioned before, the sum for (A/"^) is convergent only for integrable systems. The analysis shows that the 
shortest orbits are not providing the dominant contribution when considering the fluctuations of the particle number 
of chaotic systems. The same conclusion is reached at low temperatures for other thermodynamic functions, like 
the entropy or the specific heat. As T — > the moments of the distribution of the fluctuations of these quantities 
pick up contributions from orbits with period Tp ~ t^, whose statistical behavior is universal. Conversely, in the 
energy domain the thermodynamic fluctuations of these quantities therefore depend on the statistical properties of 
the single-particle dynamics on a scale < S. To see this explicitly in a particular case, consider the second moment of 
the entropy, given from Eq.(|9|) and Table I by 

is') = ^J 5 K(t) 42(r) . (35) 



The function K^(r) has a peak centered at r ~ of height proportional to \/T (for t/tj — > it vanishes as 
k'^(t) w —{t/tj,)^/3T). In the limit T ^ this function therefore selects in the integral (BS) arbitrarily long orbits. 
The variance of the entropy is therefore dominated by the long universal orbits whenever Tj. 3> Tmin or, equivalently, 
for temperatures kgT ^ £'c/27r^. On the contrary, for temperatures k^T « non-universal orbits in the form 

factor are selected and the result will depend on the system-specific peculiarities of the short orbits. The explicit 
dependence of (S*^) on temperature will be computed below. Similar arguments may be applied to other quantities, 
like the specific heat. We therefore conclude that those functions have, at kuT <C £'c/27r^, universal fluctuation 
distributions whose functional form depend only on the universality class of the local single-particle fluctuations on 
a scale ~ 5. 



IV. THE SECOND MOMENT IN THE tm/jv-APPROXIMATION 



There is in fact a simple way to estimate the variance in all the cases considered above, independently of the short 
or long-time origin of the dominant contributions. The approximation consists in using in Eq. ( p9| ) the appropriate 
statistical form factor (given by Eqs. (|o|) and ( |3l| ) for integrable and chaotic motion, respectively), the corresponding 



9 



value of Q'^(t) from Eq. (|3^), and to put K{t) = for r < Tmin in all cases. This is clearly an approximation, 
that we call the rmi„-approximation, since we are extrapolating the statistical behavior of the orbits down to times 
T ~ Tjnin, ignoring the short-time system-dependent structures. All the short-time structures are condensed into a 
single parameter, the period of the shortest orbit. For the fluctuations of the thermodynamic quantities dominated by 
the shortest orbits ~ whose second moment is described more accurately below, see Eq. ( ^ ) - this is clearly a rough 
approximation. On the contrary, for the entropy at low temperatures, and for the variance of the particle number in 
chaotic systems, the leading order will be well described by this approximation while the error, as we shall now see, 
is made in the correction terms. 

The virtue of the rmi„-approximation is to provide simple estimates of the size of the fluctuations, as well as of 
their dependence with chemical potential, particle number, or external parameter. It requires a minimum amount of 
information of the system, and is thus of interest in the analysis of experiments. 



A. Zero temperature. 



Setting Kj. — 1 and k'^ = in E q.(|29| ), in this approximation the results are expressed in terms of Tmin, the 
parameters a and w entering in Eq. (|33D, the shell strength q d efined in Section [I C , and the symmetry parameter 
(3 in chaotic systems. Since the integrals obtained from Eq.(P9|) are straightforward, we do not give here a detailed 
account of their computation. The results obtained are summarized in Table II. For chaotic systems with time reversal 
symmetry {(3 = 1), only the leading order terms in l/g are included. The first row of the table also describes the 
variance of the grand potential. 





Integrable 


Chaotic (3 = 1 


Chaotic /? = 2 






1 

8^ 


El( 


9 ) 






'-^) 










J^log5 + l-y 












1 

ijJ Q 

27r2 ^ 


1 


a 






1 a 













Table II: Zero-temperature second moments obtained in the T„ii„-approximation 

From these results it is possible to establish the scaling of the variances with the Fermi energy and with the system 
size L for a gas contained in a d-dimensional cavity. We assume here that the parameter x entering in the definition 
of the response is a geometrical parameter that controls the shape of the cavity (other parameters can be treated 
likewise). Since Tmin « LE^^^^, then Ec oc L^^Ep^"^. From Eqs. (|l2|) and (|l3|) (generalized to d-dimensions) we 
have (5-1 oc L'^ijf It then follows that g cx (LeI^^)'^-^. It has moreover been shown that for shape 

deformations of a cavity a oc L~^Ep^^ , and to oc L^^E^. Using these relations, we obtain the results summarized in 
Table III. The dependence with the number of fermions N in the gas has been computed assuming a constant density 
n = N/V of the gas. Therefore Ejp remains constant and, according to Eq. (p^, the size of the system increases like 
L oc A^i/'^ with the number of particles. Different dependencies are obtained when the volume is kept fixed. It is also 
interesting to compare the typical quantum fluctuations to the magnitude of the average part. For instance, since at 

constant density C7 oc TV (for any dimension) (cf (lO|), then J (U^) /U oc l/Ar(rf+3)/2<i^ 





Integrable 


Chaotic 


(c72) 


^d-3_g;(d+l)/2 ^ ^[d-Z)/d 


L-^Ep ~ N-^^d 




{LEl'y-^ ^ N'^d-l)/d 


\og{LE'J^)^logN 




j^d-3^id+3)/2 _ jy(d-3)/d 
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Table III: Zero-temperature scaling of the variances with the system size L, Fermi energy, and with the number of 
particles N in the gas in a d-dimensional cavity (keeping the density fixed). 



B. Temperature dependence of the fluctuations 



The variances depend on temperature through the function k.^. (or its derivative), contained in the function Fir) 
in Eq.(29). It would be too cumbersome to go through each function in detail. To be specific, we illustrate the 
results by computing the temperature dependence of the second moment of the entropy fluctuations, Eg. (^5|) . In the 
Tmin-approximation, Eq.(B5|) is written 



fc 



dx ( 



K{x,Xh), 



x^ Vsinhx sinh x 

where K(x,Xh) = {^I'''t)K{t = xt^), Xh = Th/tj., and K{x,Xh) = for x < Xmim where Xmin 

0.15 I 



(36) 



<S"> / kg 




0.05 



FIG. 1. Temperature dependence of the entropy fluctuations in chaotic systems, with g = 130. Full and dashed lines: Eq. ( pq ) 
(Tmin-approximation) for /3 = 1 and (3 = 2, respectively. Squares: Fermi gas in a Sinai cavity. The inset shows the behavior 
close to the origin. The long-dashed straight line corresponds to Eq.(^). 

In chaotic systems three different temperature regimes can be distinguished. 

* Low temperatures: k^T <^ 6/2tt'^. In this regime the function k'^{t) is centered at times much larger than r^. 
We therefore ignore all the structure of the form factor at finite r and take the asymptotic limit K{ 
Eq.(|6l). Then 



dx f X x'^ cosh X 
x'^ I sinh x sinh^x 



Denoting 



sinh X sinh 



the variance of the entropy in the low temperature regime grows like 



1 



((5/27r2 



(37) 



(38) 



(39) 



where 1^ « 0.153842. The variance therefore increases linearly with temperature, with a slope proportional to the av- 
erage single-particle density of states. It is interesting to remark that, up to a constant prefactor, the low-temperature 
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growth of the quantum fluctuations agrees with that of the thermal fluctuations of the gas, which are of completely 
different physical origin |l0[| . It is also interesting to compare the typical size of the quantum fluctuations to their 



average part. From Eq.(|39|) we have y (S^) — ^l^jl k g ^Jk'^TJJJJ^n^ , to be compared to the (low-temperature) 
growth of the average part S — {l/6)kg{kBT/{S/2Tr'^)). 

* Intermediate temperatures: 6/2tt'^ <C kgT <C Ec/2tt'^. In this regime the form factor may be approximated by 
its linear behavior, Eq.(^2|). Using moreover the T„iin approximation for the integral we get from Eq. (^) 



1 



{S'') = ^h kl - 



\, ( k,T 
18/3 ^ V£^c/27r2 



(40) 



where /3 = 0.141832. There is therefore a saturation of the magnitude of the entropy fluctuations in this intermediate 
regime at the "quantum of entropy" y (5*2) w /cb, with a weak negative quadratic dependence on temperature. Notice 



also the inverse dependence on the symmetry parameter /?, which is absent in Eq. (39) 



* High temperatures: kgT ^ Ec/2tt'^. Since the form factor is different from zero for times r > Tmin, in this 
limit we can approximate the function ^'^(t) by its exponential tail k'j.{t) « — (2/T)(t/tt)2 exp(— r/r^). Then, in 
the rmi„-approximation and using again the linear short-time behavior of K{t) we find 



ksT 
Ec/2n^ 



-2kBT/{Ej27r'') 



(41) 



At high temperatures compared to Ec we thus observe the expected exponential decay of the quantum fluctuations 
due to the thermal smoothing. 

The behavior of the entropy fluctuations as a function of temperature computed from Eq. ( |3^ ) for chaotic systems 
in this approximation is displayed in Fig. |l| for the two symmetry classes. 



<S > / kg' 




271 T/E 



FIG. 2. Temperature dependence of the entropy fluctuations in regular systems, with g = 1500. Long-dashed line: Eq. ([ 
(Tmm-approximation). Squares: Fermi gas in a rectangular cavity. Full line: semiclassical sum over periodic orbits, Eq.(k 
The inset shows the behavior close to the origin. 



The situation is different in integrable systems due to the different behavior of the form factor, Eq.(30). Accordingly, 
there are only two temperature regimes. 

* Low temperatures: kgT <C Ecl2iP' . In this regime the form factor is approximated by its asymptotic value Th-, 
ignoring its short-time structure. The result is identical to the low temperature regime of chaotic systems, Eq.(^9|). 
The main difference with respect to the chaotic case is that now this linear growth extends up to temperatures of 



order Ec/2t:'^ instead of 5/2tt'^. At this temperature, the typical entropy fluctuations W (5^) are of order \/gk%, which 



12 



are a factor ^ ^ 1 greater than the maximum fluctuations in chaotic systems. 

* High, temperatures: k^T 3> Ec/iir'^. As in chaotic systems, the function k'j,{t) is approximated by its exponential 
tail, and use is made of the Tmin approximation for the form factor, K{t) = for r < Tmin and K{t) = for r > Tmm- 
We get 

^^'^ = (^) e-''=-^/(^=/'^') . (42) 

Fig. 1^ shows, in this approximation, the behavior of the second moment of the entropy fluctuations in integrable 
systems. 

The second moment of the entropy fluctuations computed in the "approximation are expected to be accurate 
in the regime 2'iT^kBT ~ (5, where only long orbits (compared to Train) contribute. This is confirmed by comparing 
with numerical simulations of a Fermi gas in a cavity. Two different geometries are considered (cf Fig. In the 
first, a Fermi gas is contained inside a 2D cavity with infinite wall potential having a Sinai shape. The corresponding 
classical single-particle motion is fully chaotic. The second is similar, but the cavity has a rectangular shape. This 
guarantees an integrable single-particle dynamics. The results obtained for the numerical simulations of a Fermi gas 
in these two geometries are plotted in Figs. Q and |^ (squares). A good agreement is observed indeed in the regime 
2TT^kBT < 6 (cf the insets). 




x = r/d x = a/b 

FIG. 3. Cavities used to confined the Fermi gas. Left: Sinai (chaotic). Right: rectangle (integrable). The parameter x as 
defined in each case is varied to modify the shape (while keeping the area constant). 



As temperature increases, shorter orbits become relevant for the entropy fluctuations. For temperatures 27r^fcBT « 
Ec or higher the Tmm-approximation fails because a statistical treatment of the orbits is not appropriate any more. 
These deviations are clearly seen in Figs. and | In this range a precise description of the decay of the fluctuations 
requires to explicitly take into account the different short periodic orbits of the single-particle motion, which are 
dominant. This treatment, which requires a much m ore d etailed knowledge about the system than the simple Tmin~ 
approximation, is considered in the following section (VB). It produces in the case of the rectangular cavity the much 
more accurate description represented by the full curve in Fig. ^. 

The temperature dependence of other thermodynamic quantities may be computed analogously. For example, the 
second moment of the energy of chaotic systems shows a quadratic decrease up to temperatures fc^T « Ec/2tt'^, 
followed by an exponential tail for temperatures fc^T ^ Ecjlit^. A qualitatively similar behavior is obtained for 
integrable systems. The temperature dependence of the second moment of the particle number and of the response 
function may be computed similarly. 



V. HIGHER MOMENTS AND DISTRIBUTIONS 



A. Universal thermodynamic fluctuations 



In this case the fluctuations can be directly related to the local universal statistics of the single-particle energy 
levels. The latter correspond to an uncorrelated sequence for integrable systems, and to a random matrix sequence 
in the chaotic case. This type of fluctuations are observed in some thermodynamic functions at temperatures k^T <ti 
EJ{2t:^). 

The simplest quantity having a universal distribution at low temperatures is the particle number in chaotic systems 
with a fixed chemical potential. It has been conjectured that at zero temperature the fluctuations of the particle 
number in chaotic systems are Gaussian distributed (cf part (b) of Fig Since the variance of M was already 
computed in Section IV, the distribution is completely fixed in the range kgT -C Ec/2it'^ where the universality holds. 
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Another thermodynamic function having universal fluctuations that we consider here in detail is the entropy. When 
expressed in terms of the single-particle energy levels Ej, this function, defined by Eq.(H), takes the form 



S{x,fi,T) = fc^ 



3 



log 



T ^ 1 



_j_ ^(Ej-ti)/kE 



(43) 



At zero temperature, and as a function of the chemical potential, S{x, ^, T) vanishes everywhere except when fi = Ej, 
J = 1, 2, . . ., where it takes the value S = kg log 2. As temperature increases, peaks of width proportional to T centered 
at /i = Ej develop. At a given chemical potential /i, the entropy is therefore the superposition of the contribution 
of the different peaks. In the range kgT <^ 5/2t:'^ this superposition is only determined by the neighboring energy 
levels Ej « /x, and by their local statistical properties. We thus see, in a direct manner, that at low temperatures the 
statistical properties of the entropy depend only on the local statistics of the single-particle energy levels. The second 
moment of that distribution as a function of temperature was computed in Section IV. 

To compute the full distribution for regular and chaotic dynamics in the range fc^T <C i?c /27r^ one therefore has to 
compute the distribution of the entropy S from Eq. (^) by assuming an uncorrelated distribution and a random matrix 
distribution for the i?j's, respectively. We have done a numerical computation of those probability distributions as a 
function of temperature. The results obtained are summarized in Figs. ^ and ^. 




FIG. 4. Histogram of the entropy probability distribution for integrable single-particle motion at different temperatures. 
Dotted line: uncorrelated spectrum. Full line: Fermi gas in a rectangular box with x = (\/5 + l)/2, /i ~ (5.75 ± 0.75) x 10^, 
and g ^ 1500. (a) keT = S/{2tv^); (b) kgT = 8 x S/{2tv^); (c) kgT = 32 x 5/{2tt^); (d) ksT = 512 x 5/{2tv^) ^ {l/i)Ec/{2%^). 
In parts (a), (b) and (c) the two curves are almost indistinguishable. 

Figure |^ compares the distribution of the entropy fluctuations obtained from Eq. (^^ when the single-particle 
energy levels Ej are assumed to be uncorrelated with the probability distributions computed for a Fermi gas in a 
rectangular box. For temperatures up to 2n'^kBT « 306 « Ec/50 (part (a), (b), and (c)), the two distributions are 
almost indistinguishable and show a high sensitivity to temperature variations. For temperatures of the order Ec/2'k'^ 
or higher (part (d) of Fig. ^ the universality of the probability distribution is lost and strong deviations are observed. 
The entropy probability distribution for the uncorrelated spectrum tends to a Gaussian. 

Figure |^ shows a similar comparison but for a chaotic single-particle dynamics. It compares the probability distri- 
bution of the entropy fluctuations computed from Eq. ( ^ ) when the single-particle energy levels Ej are eigenvalues 
of a random matrix ensemble with /3 = 1, with the probability distribution computed for a Fermi gas in a Sinai 
box. Again, a remarkable agreement is found for temperatures 27r^/c£,r <C -Be, and departures are observed as E^ is 
approached. The entropy probability distribution for the random matrix spectrum tends to a Gaussian. 
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FIG. 5. Histogram of the entropy probability distribution for chaotic single-particle motion at different temperatures. Dotted 
line: random matrix spectrum (/3 = 1). Full line: Fermi gas in a Sinai box with r = 0.765 ±0.0025, fi = 900 ± 100, and g ~ 130. 
(a) fcsT = S/{2n^); (b) fcsT = 4 x S/{2n^); (c) ksT = 8x 5/{2n'^)- (d) keT = 32 x S/{2ii^) « {l/4:)Ec/{2n'^). 

For temperatures of order £'c/27r^ or higher (part (d) in Figs, ^and ||), where the universahty is lost, the moments 
of the probabihty distributions, which are now determined by the short periodic orbits, are non-universal and may 
be computed using the techniques developed below. One should keep in mind that the scale and the shell-strength 
parameter (7, set in a real system by the inverse time of flight across the system, are totally absent and therefore not 
defined in a purely uncorrelated or random matrix spectrum. 



B. Non— universal thermodynamic fluctuations 

They occur when the fluctuations are dominated by the non-universal aspects of the short time dynamics. This 
type of fluctuations are observed in all the thermodynamic quantities for temperatures k^T > _Ec/(27r^). But they 



also occur, as we saw in section [II, in cer tain thermodynamic functions at arbitrary temperatures (including T = 0). 



The Tmin-approximation used in Section IV to compute the second moment is not very accurate in this case since 
that approximation is too crude. Here we provide accurate expressions to compute all the moments of the distribution 
for non-universal fluctuations. 

Since off diagonal terms involving different short periodic orbits in products like Eq. ( |27| ) are killed by the average 
procedure, the dominant contribution to the double sum will come from the rapidly convergent diagonal terms p = p' . 
From Eq. ( p7[ ) the dominant diagonal contribution to the second moment of thermodynamic quantities with non- 
universal fluctuations, that we denote by a subscript "0", is 

{^)=2C'J2K.r- (44) 

The advantage of the Tmi„-approximation computed in the previous section compared to the more accurate expression 
Eq.(p4[) is that only information about the period of the shortest orbit is required, while Eq.(Q) involves much more 
detailed information about the periodic orbits. 
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-0.5 0.5 -0.6 0.6 -0.8 0.8 

FIG. 6. Histogram of the probability distribution of normalized thermodynamic functions for a Fermi gas in a 2D rectangular 
box: (a) U/^i^'^, (b) N/^i^'^, (c) R/^£''*'. The distributions are computed at a; = (l + ^)/2, in the window ^ = (1±0.1) x 10® 
(3 = 1500). Full line: T = 0. Dashed line: ksT = 1273 x 5/{2tt^) ~ 0.5 x £c/(27r^). 

The diagonal approximation Eq. (Q) can be generalized to compute all the moments of the oscillating quantities. 
Unlike the variance, a main difference concerning the higher moments is that the repetitions of the periodic orbits play 
now an essential role. We will explicitly go through the calculation of the third and fourth moments, the generalization 
to the higher ones being straightforward. We return to the notation of equation (|20|), keeping in mind the dependence 
on chemical potential, parameter and temperature of the variables, but not writing it explicitly. By the semiclassical 
formula the third power is written as a triple sum 

Pl,ri P2,r2 P3,r3 

where ctp-^n — fiSp-fh + Vp^^n- The product of three cosines gives 

8 cos o^p-j^ j'^i cos cXp^ cos o^p^ ^ — 2 cos (^o^p-^ -^^i I - i^o, ^ ^^P3 ■'^s ) I 2 COS {^(y,p-^ -^^i I ^^p2 -'^2 ^^P3 -'^s ) 

I 2 COS^tkp-j^ ,^1 ^~^P'2 ,^2 ^ ^~^P3 ) ^ ^ COS^tkp-j^ ,^1 ^~^P'2 ,^2 ^^p3 7^3 ) ' 

The average kills the first term in the r.h.s.. The last three terms give the same result interchanging indices. Then 

($ ) = 6 C ( ^ ^ ^ ^ ^ ^ -4pi,ri-4p2,r2-4p3,f3 C0s(Q;pj^,.j + Q!p2,r2 ~ '-'^PSii's) / ■ (46) 

\pi,ri P2,r2 P3,r3 ' fi,x 

The main dependence of this function on the chemical potential comes, as usual, from the variations of the argument 
A(j) of the cosine. In the averaging window it is enough to consider a linear approximation 

A(j) = (riTp, + r2Tp2 - rsTp^) ^ . (47) 

The wild oscillatory behavior produce by the phase variations will typically have zero average, except for those terms 
for which = 0. Assuming incommensurability of the periods of all the orbits, the only way to satisfy this condition 
is by imposing 

Tpi = Tps , r-i + rs = rg . (48) 

This is a generalization of the diagonal approximation, and illustrates the basic mechanism to calculate the higher 



moments. In the argument of the cosine in Eq. (46) only the Maslov indices survive because the equality of the 
periods implies the equality of the actions. The above constraints simplify the multiple sums and finally a compact 
expression for the third moment is obtained 



oc oo 



(^o) — 6 C ^ ^ ^ ^ ^ ^ •^p,ri •^p.r2 •^p,ri-\-r2 ^^^{,^p^ri ~t~ ^p^r2 ^p^Ti-\-r2^ ■ (4^) 



p ri = lr2 = l 
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When the fluctuations are dominated by the short orbits, the convergence of this series is guaranteed by the same 
arguments that showed the convergence of the second moment. The sum ( ^9|) is in general different from zero. This 
produces a finite skew giving a non-symmetric probabihty distribution. 



k 










T = 




T = 


Ti 7^0 


T = 


Ti 7^0 


2 


(1.81 ± 0.01) X 10-2 
1.83290 X 10-2 


(1.55 ±0.01) X 10-2 
1.57853 X 10-2 


(2.83 ± 0.05) X 10-2 
2.8477 X 10-2 


(1.63 ±0.03) X 10-2 
1.65251 X 10-2 


(4.25 ±0.10) X 10-2 
4.3148 X 10-2 


(2.84 ±0.06) X 10-2 
2.89043 X 10-2 


3 


(5.0 ±0.4) X 10""' 
5.31694 X 10"'' 


(3.5 ± 0.5) X 10-"' 
3.81036 X 10-" 


-(1.3 ±0.2) X 10-3 
-1.35453 X 10-3 


-(5.5 ± 0.7) X 10-"' 
-5.85974 X 10-"' 


(2.4 ±0.3) X 10-3 
2.44332 X 10-3 


(1.6 ±0.1) X 10-3 
1.58641 X 10-3 


4 


(7.1 ±0.2) X 10""' 
7.3133 X 10""' 


(4.95 ±0.05) X 10-"* 
5.13214 X 10-"* 


(2.3 ±0.1) X 10-3 
2.33807 X 10-3 


(6.8 ± 0.3) X 10-"' 
7.14034 X 10-"' 


(5.0 ±0.3) X 10-3 
5.16324 X 10-3 


(1.9 ± 1.0) X 10-3 
1.99110 X 10-3 



Table IV: Moments of the distribution of some thermodynamic functions with non-universal fluctuations for a Fermi eas in a 
rectangular cavity. In each row, the upper values are the numerical results obtained from the distributions shown in Fig. rJ, while the 
lower values are obtained from the semiclassical sums over periodic orbits. The temperature Ti = ksT = 1273 X 5/{2n^) is the same as 

in Fig. |6| 

The fourth moment is calculated in a similar way. The product of four cosines gives eight terms. One of them, 
containing the sum of all the actions, has zero average. The other terms, arranged after interchanging indices, give 

Pl,ri Pi,ri 

+ 3 COs(api,ri +ap2,r2 - Otp3,r3 - OLp^.Ti) ] ^ ■ (50) 

We restrict again to those terms having exactly zero phase variation, neglecting quasi cancellations of long orbits. Then 
for the first cosine in the r.h.s. this is equivalent to tit^-^ + r2Tp^ + r^Tp^ — T^Tp^ = 0, and by the incommensurability 
of the periods the only solution is 

Pi = P2 = P3 = P4 , + r2 + ra = r4 . (51) 

As for the third moment, this condition gives a simple sum over primitive periodic orbits. The three indices ri,r2 
and r3 are free, and r4 is determined by them. 

For the second cosine the corresponding equation riTp^ + r2Tp^ — r^Tp^ — r^Tp^ = admits two different solutions. 
The first one is 

Pi = P2 = P3 = P4 , 7-1 + r2 = ra + r4 . (52) 

The difference is that only ri and r2 are free. Because the repetitions are positive integer numbers, varies between 
1 and ri +r2 — 1, and is determined by the constraint. The other solution, which has a weight two because of the 
possibility of interchanging indices, is 

P3 = Pi , P4 = P2 , r3 = ri , ri=r2 . (53) 

This reduces the expressions to double sums over two different primitive periodic orbits, pi andp2, and the repetitions 
ri, r2. To eliminate the inequality between pi and p2, we can add and subtract a term where the two indices are equal. 
With this term the double sum, both in primitive periodic orbits and repetitions, factorizes giving a term directly 
proportional to the square of the second moment, in its diagonal approximation. Finally we arrive at a convergent 
formula for the fourth moment 

oc 

($q) 2 C ^ ^ 4 ^ ^ -^PjTi •^p,r2 -^p.rs -^p,ri -t-r2+r'3 ^^^{^p,ri ~^ ^p,r2 ~^ ^p^r^ ^p,ri+r2+r3^ 
P ri,r2,r3 = l 

oc ri+r2 — 1 

± 3 ^ ^ ^ ^ -^p.ri •^p,r2 -^p,r3 •^p,ri4-r2 — rs ^^^{,^p.ri ~^ ^p,r2 ^P^T3 ^P-T\+r2—r3^ 
ri,r2 = l r3 = l 
oo 

- 6 K,r, K,r^ +3 {^Y ■ (54) 

'ri,r2 = l 
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The last term 3($^)^ is the result for the fourth moment if the distribution would have been Gaussian. The remaining 
terms in Eq. (|5J) produce deviations from that behavior (the excess of the distribution). 

The generalization to higher moments is straightforward. The formulae are more complex, but they are always 
written as a sum over one primitive periodic orbit index, the repetitions originating from different constraints, plus 
terms coming from the previous moments (e.g. for the fifth moment there is a term proportional to(<i>^) ($"^)). 

Figures ^ and illustrate these results. They show the distribution of the quantum fluctuations of the total energy, 
particle number and response of a Fermi gas contained in a 2D rectangular and Sinai cavity, respectively. The 
fluctuating part of the thermodynamic functions is normalized by an appropriate power of the chemical potential. 
This normalization eliminates their zero-temperature energy dependence (cf Table III with d = 2). The probability 
distributions are normalized by setting their area to one. The parameter x used to compute the response controls 
the shape of the cavity (cf Fig. H), while the area is kept constant. R is therefore a "quantum" force, related to the 
variations of the energy of the gas produced by the quantum shell corrections. 

The first four moments of the probability distributions shown in Fig. ^ have been computed. The numerical results 
obtained, at the two different temperatures, are compared in Table IV to the first four moments calculated using the 
analytic expressions obtained above based on the periodic orbits. The overall precision obtained is of the order of 2% 
(in the worst case), and is within the numerical uncertainties. For comparison, the rescaled second moments of the 
energy, particle number and response computed at T = using the Tmm-approximation are 1.74 x 10^^, 3.2 x 10^^ 
and 3.98 x 10~^, respectively. They have an accuracy of 4%, 13% and 6% with respect to the numerical values of 
Table IV, respectively. 




FIG. 7. Histogram of the probability distribution of normalized thermodynamic functions for a Fermi gas in a 2D Sinai box: 
(a) f//^i^/*, (b) Mlii^l^, (c) Rl[£'l'^. The distributions are computed in the windows r = 0.765 ± 0.0006, ^ = 750 ± 250. FuU 
line: T = 0. Dashed line: ksT = 0.2 x Ec/{2n'^). 

We have also computed from Eq. (^) the temperature dependence of the variance of the entropy fluctuations for a 
Fermi gas in the rectangular box using the shortest orbits. The result is represented by the full line in Fig. ^ which 
shows a very good agreement with respect to the numerical values computed directly from the Fermi gas. Such a good 
agreement is found even at temperatures of order i5 because, in integrable systems, the "diagonal" approximation used 
above is almost exact. This is not expected to occur in chaotic systems. We were not able to make similar comparisons 
for the chaotic Sinai cavity because we haven't computed the periodic orbits for that system. 

The fluctuations of the particle number in the chaotic case. Fig. 0(b), shows at T = the universal Gaussian 
behavior discussed in section V.A. At higher temperatures the universality is lost. 

VI. CORRELATION FUNCTIONS 

When the parameter x or chemical potential /i are varied, the autocorrelation function of a thermodynamic quantity 
$ is defined as 

C$(2:,Ai) = ($(a;o-x/2,Mo-/^/2,T) $(xo + 2:/2,^o + /V2,r)\ . (55) 
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In the following, we assume that the regular or chaotic nature of the single-particle motion in the gas is unchanged 
when the parameters (xq, fio) vary. When $ is expressed as a sum over the periodic orbits using Eq.(pO[), C$ takes the 
form of a double sum over periodic orbits labelled by the indices p and p' . As usual, only the term involving one half 
the cosine of the difference of the actions gives a non-zero average. The amplitudes Ap and actions Sp are evaluated 
at (fiQ — Xq — x/2), while those with index p' at (/xg + /^/Z, + a;/2). In the semiclassical regime these functions 
will have a slow variation in a quantum scale, and therefore to leading order the parameter dependence can be ignore 
in the prefactors, and only a linear expansion of the actions is taken 

Sp{^lQ ± ^, xo ± |) w Sp{^I.o,xo) ± ^ ± ■ 

The autocorrelation is now written 



oc 



C^ix,fi)=2C' 2^ Ap,rAp<,r' COS 

\p,p' r,r' — l 



rSp - r'Sp' [rTp + r'Tp>) {rQp + r'Qp,) 



M H t;^ X 



2h 2h 



(56) 



where all the classical functions are evaluated at (a;o,/io). The average over a small window Aa; around xq restricts 
the sums to orbits having approximately the same Qp, while the energy average imposes the same condition over the 
periods. 

To proceed further there are two main paths, corresponding to universal and non-universal fluctuations. For the 
latter, since short orbits dominate, the diagonal approximation furnishes the main contribution to the autocorrelation, 

C*,o(x,^) = 2C2 Y^f^^l^^osf'-^ + '-^] . (57) 

p r—l ^ 

In this equation the shortest periodic orbits have the largest prefactors Ap.r and the lowest frequencies inside the 
cosine, giving rise in general to an erratic and non-decaying oscillation as a function of the parameters. The typical 
period of these oscillations is dominated by the shortest periodic orbit. The correlation lengths, both in chemical 
potential and external parameter x are estimated making the phase variation equal to 27r giving 

= (58) 
'min I ^min \ 

In the presence of universal fluctuations, the treatment of Eq.(p6|) is different. Since long orbits give an impor- 
tant contribution, the diagonal approximation is not sufficient and a more accurate evaluation is needed. General 
expressions for the autocorrelation of the gas were given in Ref. Ordering the orbits by their period, and taking 
into account the restrictions imposed by the averaging procedure, the autocorrelation of a thermodynamic quantity 
$ takes the form 

C*(a;,M) = Y^ ^ cosifiT/n) K{r,x) F{t) . (59) 

The constant C and the function F{t) are defined as in Eq.(p9|), to which Eq.(p9|) reduces when x — ^ — {K{t, x — 
0) = K{t)). The function K{t, x) is a generalization of the form factor: it is the Fourier transform of the parametric 
two-point correlation function (p(/io — ij,/2,xq — x/2)p{po + ^/2,xo + x/2))^g with respect to ii. Semiclassically its 
expression is given by Eq.(^) with the replacement cos[{Sp — Sp>) /h] cos[{Sp — Spi) /% + [Qp + Qpi) /2K\. For 
chaotic systems and for times r ^ Tmim K{t, x) has been conjectured to be a universal fimction depending only on 
the symmetry parameter /3 and described by random matrix theory p^ . 

In many cases this integral may lead to results which are determined only by the universal, long time aspects 
of the dynamics. Consider for example the autocorrelation of the entropy. Then F{t) = k^^(t) in Eq.(p9|). At 
temperatures k^T <^ Ec/2t:'^ the function k^^(t) selects only times much longer than Tmin- Then the integral in ( [59| ) 
is insensitive to the system-dependent features of K{t,x) present at short times, and will therefore take a universal 
form depending only on the temperature, the mean level spacing and the typical slope of the single-particle energy 
levels at Fermi energy p4|j20[ | . For (fully) chaotic systems we therefore need the random matrix theory result for this 
function. Changing for simplicity to the scaled variable t = t/th, using the definition K{t, x) = {1/th)K{t — Th t, x) 
and scaling the parameter x according to the prescription given in pO[ , for chaotic systems without time reversal 
symmetry for instance (/3 — 2) we can Fourier transform with respect to energy the parametric two-point correlation 
function obtained in |24]. We get 
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Krmtd, X) I ^^2^.^.^. s\nh{2^^xH)/2^^xH t>l. ^ ^ 

Inserting this in Eq.(|5^) the autocorrelation of the entropy may therefore be computed for temperatures well below 
Ec/2t:''. 



VII. DISCUSSION 



It follows from Eqs.(|12[) and (13) that the three "intrinsic" energy scales that characterize the Fermi gas, namely 
Ep, Ec and 5, are determined by the number of fermions in the gas (taking into account degeneracies due to, e.g., 
spin), by their mass, and by the volume of the cavity (or, alternatively, /x, m, and V fix N , Ec and d). The same 
is true for the adimensional shell-strength parameter g = Ec/6. They are independent of the precise shape of the 
system, and are therefore independent of the regular or chaotic nature of the single-particle dynamics. 





Metal particles 


Nuclei 


N 


20 - 5000 


25 - 250 


r = roN'^/'^ 


0.5 - 3.5 nm 


3.2 - 7 fm 


n = NjV 


30 electrons nm^^ 


0.17 nucleons fm""^ 


kpT 


10 - 67 


11 - 25 


Ep 


35000 K 


37 MeV 


E ''^^ 
kp-r 


22000 - 3500 K 


26.5 - 12.3 MeV 


^^2Ep 
3/V 


1200 - 5 K 


1 - 0.1 MeV 


37riV 
^ ~ 2kpr 


20 - 760 


27 - 125 



Table V: Main features of metal particles and atomic nuclei as a function of particle number, 
ranging from 20 to 5000 and from 25 to 250, respectively. 

It is useful to recall their typical magnitude as well as other relevant physical parameters for two important systems, 
metallic particles and atomic nuclei. Table V summarizes the main features. The number of particles is N, electrons 
or nucleons respectively. It typically ranges from 20 to several thousands in metal particles, and from 25 to 250 
in nuclei. For definiteness, we take the window 20-5000 for metallic particles (although values up to 20000 can be 
reached experimentally p^). The values indicated in the table are for these two windows, respectively. The radius 
r = rpiV^/^ fixes the system size as a function of N. For metal particles vq (which is sometimes denoted r^) ranges 
from 1 to 3 A. In the table we have used the typical value Ts = 2 A. For nuclei tq ~ 1.1 fm. 

Other relevant experimental systems, with much larger typical sizes, are the two-dimensional (2D) quantum dots. 
For instance, the electron gases contained in microscopic 2D cavities created in high mobility GaAs heterojunctions 
have a typical electronic density of ti 5 x 10^^ cm^^ and Fermi energy Ep — 18 meV (see, e.g., Ref. [^). For 
typical quantum-dot sizes L « 0.1 - 5 fj,m, the other characteristic energy scales vary in the range Ec ~ 6.4 - 0.128 
meV and S ~ 0.36 - 1.44 x 10~^ meV. This gives an effective number of fermions in the last shell in the range g = 18 
- 890. 

The table illustrates the variations in the energy scales with the number of particles. Ec is comparable to Ep for 
small grains and nuclei. As the system grows, the three energy scales become however well separated, thus improving 
the applicability of the semiclassical methods. 

As we have shown in the previous sections, the parameter g which, according to Table V, may take quite large 
values in real systems, plays an important role in setting the scale of the fluctuations of thermodynamic functions. 
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As summarized in Table II, it determines the relative value of the fluctuations in integrable systems as compared to 
chaotic ones for quantities like the energy or the response of the gas. This amplification factor is due to the special 
organization of periodic orbits in regular systems, which form one-parameter families, to be compared to the isolated 
character of periodic orbits for chaotic motion. Among the chaotic systems the variance is twice bigger in systems 
with time reversal symmetry. 

It is instructive to compare the fluctuations of the energy and of the response function of the whole Fermi gas with 
the corresponding fluctuations of a single-particle energy level located in the neighborhood of the Fermi energy. When 
an external parameter is varied, the fluctuations of the energy of a single-particle level located at E ^ Ep may be 
computed from the fluctuations of the Fermi energy of a system with a fixed number of particles N . The number of 
particles is related to the Fermi energy through the spectral counting function, N — M{Ep). Writing Ep ~ Ep + Ep 
and also splitting M into its smooth and oscillatory parts, expanding when appropriate by assuming Ep <^ Ep, and 
ignoring second order terms in the fluctuations we arrive at = M{Ef) + p{Ep)Ep +N'{Ep). But Ep is precisely 
defined in order to satisfy Af(Ep) = N. The fluctuations of a single-particle level, estimated by Egp ~ Ep, are 
therefore 

E,p = -M(Ep) S (61) 

with the variance 

{E%) - (A^2) ■ (62) 

This relation allows, together with the variance of the particle number summarized in Table II, to connect the energy 
fluctuations of the gas with those of the single particle levels. To leading order in g we get 

^ Td'^iEsp) integrable , 



127r2- 

iU') = < ^ (63) 



2 



1 9 
87r2 logg 



{E^ ) chaotic 



To establish a similar connexion for the response function we must first determine the local fiuctuations of the 
response of the single-particle energy levels 

Rsp^'-§, (64) 

where Ej « Ep. It was shown in Ref. that in chaotic systems the variance of Rgp (with respect to a mean which 
is subtracted) is 

/52 \ 2a 2q; 1 

(Rsp) = -0^- ' (65) 

f-J I H y ' rain y 

where a was defined in Eq.(|3^). A similar computation for integrable systems gives, using the corresponding relation 
in Eq.(|3|) 

{Rip) = ^- (66) 
Using these relations together with the results for the variance of the response of the gas (i?^) from Table II we find 

{R^) = ^9{R%) (67) 

for both regular and chaotic dynamics. 

The global picture emerging from these relations is quite instructive. Contrary to naive expectations, Eq.(|6^) shows 
that the typical fiuctuations of a single-particle level located at Ep are not of order 6, but are "amplified" by the 
fluctuations of the particle number. Since the latter are larger in integrable systems compared to chaotic ones (by a 



relative factor y^g/logg ), then y {Esp)int/{Egp)ch oc y/g/\ogg. Eq.(|63|) shows, moreover, that the fiuctuations of 
the energy of the gas are, in turn, much larger than the local variations. Similar considerations hold for the response. 



We see from Eq.(|67|) that the typical response of the gas, Rtyp = y (-R^), is ^/g times greater than that of the single- 
particle levels located at Ep, irrespective of the integrable or chaotic nature of the dynamics. One should however keep 
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in mind that (R'^p) is bigger in integrable systems with respect to chaotic ones by a factor of order g, and therefore 
the overall fluctuations of the response of the gas are more important if the dynamics is integrable. When computed 
in the particular case of the persistent currents in mesoscopic rings, I — —c d^l/dcj), where the external parameter x is 
now the magnetic flux threading the ring, the formula (i?^) — a/ (n'^Tmin) in Table II valid for chaotic systems gives 

for the typical current Ityp = \J {P) = {■s/2c/t:)EcI4>{) (where c is the speed of light and (po the flux quantum), in 
agreement with the result obtained in [2^ using a Green's-^function approach, and with the y/g amplification obtained 
for integrable systems (see also p^). 

All this results clearly show the close connexions existing between the fluctuations of the gas and the parameter g. 
Remember that g = Ed 5 is the number of fermions contained in the last shell, i.e. in a window of size Ec below the 
Fermi energy. This number of particles has to be compared with the total number of fermions N = {2/d)Ep/5 (in 
d dimensions). Eqs. (^^ and (^^ show that neither the total number of particles, nor an individual single-particle 
level (located at Ep), are responsible for the fluctuations of the gas, but rather the intermediate number g, with 
typically 1 <C 5 ^ iV. The quantum fluctuations of the gas therefore originate from the contributions of the g 
particles located in the last shell, as already pointed out in Ref. However, the way these particles contribute to 
the fluctuations depend on the quantity considered. Eq. (m) shows that they contribute "coherently" in the case of 



the total ground-state energy, since Utyp — y (U"^) is proportional to g (ignoring the logg denominator in the chaotic 
case). In contrast, their contributions to the response add up "incoherently" (as in a random walk), as shown by 
Eq. (|67|). Physically, this may be interpreted as the independence of the slopes dEj/dx of the last g single-particle 
levels. 



The situation is more subtle for quantities like the entropy. As was shown in section EVB, the temperature 
dependence of the variance for ksT <C S/2tt'^ is unique, independent of the nature of the single-particle dynamics 



(and given by Eq. (|39|)). However, at higher temperatures differences are observed. The maximum of the typical 



value of the entropy fluctuations, Styp — y (S'^), is of order kg for chaotic systems, and it is reached at 2'K^kBT w 5, 
whereas for integrable systems Styp grows up to values oc gk^. The maximum value in integrable systems is reached 
at temperatures 2t{^ kgT m Ec- 

There are several distinct applications of the present results. Among them, we mentioned the study of the electronic 
contribution to the mechanical force in experiments where metallic nanocontacts are pulled [29|] , and the fluctuations 
of the energy of nuclei as a function of the number of nucleons. Progress made in both directions will be published 
elsewhere. In a related work, the distribution of the energy of a Fermi gas whose single-particle energy levels are 
given by the imaginary part of the complex zeros of the Riemann zeta function was recently investigated pO| . This 
fictitious fermionic system was named "the Riemannium". Aside its mathematical interest, it serves as an explicit and 
important test of our results relating the probability distribution of thermodynamic functions to periodic orbit theory. 
This is because there exist strong connexions between the quantum theory of chaotic systems and the complex zeros 
of the Riemann zeta function interpreted as spectral eigenvalues. A high precision agreement was found between the 
numerical results and the moments of the probability distribution of f2 computed from the results of Section V A for 
the Riemannium. 



VIII. CONCLUDING REMARKS 



Our results illustrate the strong connexions and deep relationships that exist between, on the one hand, the 
thermodynamic quantum fluctuations of confined Fermi gases and, on the other, the statistical theories of level 
fluctuations and periodic orbit theory. 

The regular or chaotic nature of the single-particle motion imprints the single-particle spectrum in two different 
ways. The first one occurs at the scale of the mean level spacing 5, and produces different (but universal) statistical 
fluctuations. The second acts on a much larger scale E^^ set by the inverse of the time of flight across the system. A 
bunching of the single-particle levels produced by the short periodic orbits is observed on this scale, whose intensity 
depends on the regular or chaotic nature of the motion. 

The main theme of our investigations has been to determine how these two features of the single-particle spectrum 
influence the probability distribution of the quantum fluctuations of the different thermodynamic functions of the 
gas, and to establish their temperature dependence. We find a rich variety of phenomena and of different regimes. In 
some cases the thermodynamic quantum fluctuations directly reflect the universality of the single-particle spectrum 
on the scale 5. Their corresponding probability distributions are universal functions, whose shape depends on the 
quantity and temperature considered, and on the regular or chaotic nature of the dynamics, but not on the precise 
shape or other specific properties of the system. In other cases, the fluctuations are, in contrast, totally determined 
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by the long-range modulations of the single-particle spectrum on scales Ec- They are therefore insensitive to the 
universality present on the scale 5. For some thermodynamic quantities, this type of fluctuations are dominant even 
at zero temperature. The shape of the distributions is, in this case, system-dependent and therefore non-universal. 
The moments may be computed from the short periodic orbits of the single-particle dynamics. At temperatures of 
order Ec or higher, all the thermodynamic fluctuations, even those dominated at low temperatures by the universal 
features of the single-particle spectrum, fall in this second class. In all cases and for the different regimes, we have 
shown how the probability distributions can be explicitly computed in each case. 

The variance of the thermodynamic functions has also been investigated using a simple approximation, that takes 
only into account the regular or chaotic nature of the single-particle dynamics and the time of flight across the 
system. It has the advantage of providing a good qualitative description for the different thermodynamic quantities 
as a function of temperature, number of particles, etc, while requiring a minimum amount of information about the 
system. Although more accurate descriptions can be made if more information is available, this approximation is 
clearly of interest in experiments where the parameters controlling the system, like its shape, are not well known. 

A closely related problem is the motion of quasiparticles in potentials with bulk disorder. This case has been 
extensively treated in condensed matter physics in particular as a model for the electronic properties of metals and 
semiconductors. Compared to the ballistic motion considered here, the disorder introduces a new scale, the elastic 
mean free path £e between impurities. When ig is much smaller than the system size L (and L is much smaller than 
the localization length), the single-particle motion across the system is diffusive. The corresponding energy scale is 
the Thouless energy, E^ — hD/L^, with D the diffusion coefhcient. Similarly to the influence of the (short) periodic 
orbits in the ballistic case, the diffusive motion produces long-range correlations in the single-particle spectrum 
However, the nature of the correlations are different in both cases. Therefore, the statistics of the quantum 
thermodynamic fluctuations that are dominated by the long range correlations are going to be different. An analogous 
study for diffusive systems to the one made here in the ballistic case is therefore necessary. In contrast, it has been 
shown that the local fluctuations of the single-particle levels on a scale 6 are, in the metallic regime of disordered 
systems, also described by random matrix theory |^^. Therefore, our results concerning the quantum fluctuations of 
thermodynamic quantities in the universal regime apply also to disordered systems. 

The quantum fluctuations considered may be compared to other fluctuations, the thermal ones. Although the latter 
are inherent to any thcrmodynamical treatment, their physical origin and their experimental manifestation is very 
different from that of the quantum fluctuations here considered. 

Many experiments have verified the validity and accuracy of the physical picture obtained from mean-field, single- 
particle approximations to many-body systems. However the nature of the modifications induced by the (residual) 
interactions in the probability distribution of the thermodynamic fluctuations remains an open problem. 

Acknowledgements: we are in debt with C. Schmit who kindly bring us the numerical spectrum of the Sinai 
cavity. 
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